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The universal quantum computer]!] is a device capable of simulating any physical system[2] and 
represents a major goal for the field of quantum information science. Algorithms performed on such 
a device are predicted to offer significant gains for some important computational tasks 3 . In the 
context of quantum information, "universal" refers to the ability to perform arbitrary unitary trans- 
formations in the system's computational space [2J. The combination of arbitrary single-quantum-bit 
(qubit) gates with an entangling two-qubit gate is a gate set capable of achieving universal control of 
any number of qubits 5, 6, 7 , provided that these gates can be performed repeatedly and between 
arbitrary pairs of qubits. Although gate sets have been demonstrated in several technologies [Sj, 
they have as yet been tailored toward specific tasks, forming a small subset of all unitary oper- 
ators. Here we demonstrate a programmable quantum processor that realises arbitrary unitary 
transformations on two qubits, which are stored in trapped atomic ions. Using quantum state and 
process tomography [9], we characterise the fidelity of our implementation for 160 randomly chosen 
operations. This universal control is equivalent to simulating any pairwise interaction between spin- 
1/2 systems. A programmable multi-qubit register could form a core component of a large-scale 
quantum processor, and the methods used here are suitable for such a device [TO]. 



Computers are useful because they are versatile. 
Changing the problem to be solved amounts to reconfig- 
uring inputs to the processor, that is, to reprogramming 
it. In a classical computer, a program is ultimately de- 
composed into sequences of operations implemented with 
logic gates. The explosion of interest in quantum infor- 
mation science coincided with the realisation that a sim- 
ilar decomposition exists for quantum processors [H[5J|5]; 
arbitrary operations on a multi-qubit system can be bro- 
ken down into sequences of discrete operators - "quan- 
tum gates". As with its classical counterpart, a pro- 
grammable quantum computer is more versatile than one 
designed for a fixed task. 

Ease of implementation can favor certain decomposi- 
tions of quantum operations, for example, those based 
on arbitrary single-qubit gates and a single entangling 
two-qubit gate[5l |6]. Since realising gates acting on 
two or more qubits tends to be more experimentally 
challenging [5], much attention has been focused on us- 
ing them optimally in the creation of entanglement [TT] 
and on finding decompositions minimizing the number 
of times they are applied H21 US HH [15] . Some well- 
chosen operations can be performed with two or fewer ap- 
plications of two-qubit gates, but these form an infinitesi- 
mal subset in the space of two-qubit operations |13j. For a 
two-qubit system with a maximally entangling two-qubit 
gate, three applications of this gate, when combined with 
additional applications of arbitrary single-qubit gates, 
are sufficient for universality [12 . Here we present and 
characterise a universal quantum processor that can pro- 
duce any desired two-qubit unitary transformation when 
programmed with 15 classical inputs [TT1 113] . 

The decomposition of a given operation depends on 
the available gate sets. Our choice of a universal gate 
library consists of single-qubit gates and one maximally 
entangling two-qubit gate. The single qubit gates are 
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FIG. 1: Components for universal computation. Cir- 
cuit diagrams for arbitrary unitary transformations on a, one 
and b, two qubits. The operations for each circuit are im- 
plemented from left to right with each line representing one 
qubit. Part a also indicates the decomposition employed for 
R(6, (/>). The dashed box in part b contains the three degrees 
of freedom a, f3, 5 that determine the two-qubit operation's 
local equivalence class. The brackets highlight the decompo- 
sition of the two-qubit operation U as described in the text: 
U = (C® D) ■ V ■ {A® B). 



rotations 

R(8, (f>) = exp[— iO (cos <j> a x + sin cj) o- y ) /2] (1) 
RM = exp(-i<p z a z /2) (2) 

in the computational basis |1) = (1, 0), |0) = (0, 1). Here, 
cr x ,(jy,a z are the Pauli matrices. The single-qubit gates 
have variables 9, <f>, <f> z that can take any value from to 
2ir. The entangling two-qubit gate is 



G = e-"/ 4 exp( ~a z ®a z ) (3) 



and operates on the basis |11) , j 10) , 1 01) , |00). Here, 
indicates the tensor product, and \ij) = \i) ® \ j). 
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With this gate library, the circuits shown in Fig- 
ure 1 can be used to implement arbitrary unitary trans- 
formations on (a) one qubit and (b) two qubits. The 
single-qubit operation is characterized by three degrees of 
freedom and may be decomposed into the matrix product 
R z {4>z) • R{6,4>). Each two-qubit unitary transformation 
is described by 15 degrees of freedom [TT1 113] . The group 
of all such transformations can be divided into subsets 
that are equivalent up to single-qubit operations. Such 
subsets are called local equivalence classes jH] because lo- 
cal operations can map among all members of the class. 
Each local equivalence class can be described by three 
parameters 11J. Given a unitary transformation U, we 
decompose it into U = (C <£> D) ■ V ■ {A (g> B). Here, V 
is in the same local equivalence class as U and is in a 
special form that requires fewer gates on our processor. 
A and C are single-qubit operations on one of the qubits, 
and B and D are single-qubit operations on the other. 

Determination of the 15 single-qubit-gate parame- 
ters is facilitated by working in the so-called "magic" 
basis [TT]: 

-1(|11> + |00>),-^(|11)-|00>), 

V2 V2 (4) 

~(|10> + |01>),~(|10)-|01>). 

This basis amounts to the Bell states with specific global 
phases, and we take advantage of two of its convenient 
mathematical properties. First, single-qubit operations 
with unit determinant are given in this basis by real ma- 
trices that are orthogonal [16] . Second, two-qubit oper- 
ations u and v in SU(4) are in the same local equiva- 
lence class if and only if uu T and vv T have the same 
eigenvalues [16 (it T is the matrix transpose of u, and we 
use lower-case letters to represent matrices in the magic 
basis). The decomposition of a given operation U follows 
a two-step procedure analogous to that in Ref. [15] . (See 
Methods for details and Supplementary Information for 
examples.) Briefly, we first match matrix eigenvalues to 
find a special clement V in the same local equivalence 
class as U (i.e. find a, (3, S, the three degrees of freedom 
in Figure lb's dashed box). Second, we manipulate real, 
orthonormal matrix eigenvectors to find the four remain- 
ing single-qubit operations required to map between V 
and U. Because the magic basis properties rely on unit 
matrix determinants, we can implement operations only 
up to a global phase. Global phases exactly vanish in 
any observable quantity, so this restriction has no physi- 
cal relevance. 

We implement the quantum circuit with trapped 
ions using techniques applicable for scaling to a larger 
svstem[10] . Each qubit is stored in a pair of en- 
ergy eigenstates in the 9 Be + 2s 2 S , 1 / 2 hyperfine mani- 
fold. The qubit basis states can be transferred be- 
tween different pairs of the eight hyperfine levels [TO] . 
The qubit spends the majority of its time stored in 



the (|1),|0)) = (\F = l,m F = 0), \F = 2,m F = l)) 
"magnetic-field-independent" manifold, for which the en- 
ergy splitting has zero first-order dependence on the mag- 
netic field at our chosen value of 0.011964 T, leading to 
long coherence times (15 s has been measured for a pair 
of states with similar second-order field dependence [T7]). 

The two 9 Be + ions are stored simultaneously with two 
24 Mg + ions in a six-zone linear Paul trap [101 ESI US]; the 
ions form a linear chain along the axis of weakest confine- 
ment. Coulomb repulsion couples the motion of all four 
ions such that laser cooling the 24 Mg + ions sympatheti- 
cally cools all of the ions without destroying the quan- 
tum information |10| [19] stored in the 9 Be + ions. The 
collective motion of the four ions can be described as the 
sum of 12 normal modes, four along each of the principal 
axes. The two-qubit gate uses two modes involving mo- 
tion along the axis of the ion chainJTO] [19] . We spectrally 
address the 24 Mg + ions to Doppler and resolved-sideband 
cool these modes to near the quantum ground state of 
motion [TU]. Since the ion spatial order affects the mode 
frequencies, and because both resolved-sideband cooling 
and the two-qubit gate require spectral addressing of the 
modes, we deterministically initialize the ion order [19] 
to 9 Be + - 24 Mg + - 24 Mg + - 9 Be+ at the beginning of each 
experimental sequence. 

State preparation and measurement are performed us- 
ing resonant laser light that couples the 9 Be + 2s 2 S , 1 / 2 
states to the 2p 2 P 1 / 2 and 2p 2 P 3 / 2 statesj^O]. Projec- 
tive measurements in the single-qubit computational ba- 
sis |1) , |0) utilize a cycling transition [TU] HPJ. Measure- 
ments in other bases are made by first rotating their 
eigenvectors into the computational basis. A single de- 
tection apparatus sequentially measures the two qubits 
independently [TO] . 

The universal gate set above is implemented with laser- 
induced stimulated- Raman transitions [20] . The two- 
qubit gate G is a geometric phase gate [TO] EJJ. The 
single-qubit gate R(0, 4>) can be produced by driving res- 
onant Rabi oscillations between the qubit states, where 
the angle 6 is controlled by the laser pulse intensity and 
duration. The phase <j> is set by the phase difference 
between the two Raman light fields at the ion relative 
to the qubit phase [50]. It is controlled via the phase of 
a radio-frequency (RF) potential applied to an acousto- 
optic modulator (AOM). The single-qubit gate R z ((j) z ) 
advances the qubit phase by (f> z relative to that of the Ra- 
man light fields. It is implemented by retarding the phase 
of an AOM's RF for subsequent laser pulses. In order to 
individually apply R(0, <fi) to each qubit, time-dependent 
electric potentials [22 divide the four-ion chain into two 
9 Be + - 24 Mg + pairs and transport them to zones sepa- 
rated by 240 /im. The four applications of R(0, </>) per 
qubit in Figure lb require four such ion separations and 
subsequent recombinations - a total information trans- 
port of nearly 2 mm per qubit. An AOM placed before 
the trap directs the laser beams to the relevant ion. 



FIG. 2: Diversity of states. Examples of the density matrices created by applying four of the randomly chosen unitary 
transformations to simple product states. The bar heights and colors show the absolute value of the output density matrices, 
with the experimental data shown solid and the theory shown semi-transparent. These states have fidelities / = 0.81, 0.80, 
0.75, 0.80, relative to the ideal states. The Supplementary Information gives the input states and operations used in these 
examples. 



The requirement that single-qubit-gate inputs take any 
value from to 2ir prompts a further decomposition of 
R(8,4>). The RF potentials that control 4> an d 4>z are 
generated by a stable oscillator whose phase is easily con- 
trollable. The variable 8, however, depends on the laser 
pulse's intensity, which for technical reasons is not con- 
stant for the duration of the pulse. Rather than calibrate 
this for arbitrary 8, we calibrate a single value, 8 = tt/2, 
and decompose R(8, cf>) into 

R{8,ct>) = |) • RM ■ R^<t>-\) • (5) 

In this way, all 15 inputs to our universal circuit are 
controlled by shifting the phase of a control oscillator 
relative to the qubit. The number, duration, and spacing 
of the laser pulses are identical for every U. 

To demonstrate the ability of the processor to generate 
arbitrary unitary transformations, we program it with 
160 different randomly chosen operations distributed in 
SU(4) according to the Haar measure[23]. (The proba- 
bility distribution given by the Haar measure is a uni- 
form distribution in the space of unitary matrices.) To 
characterise our implementation of the 160 operations, 
we apply each to one of 16 input states formed by the 
tensor products of |0>, |1), |+> = + |0)) /a/2, and 
= — i |0)) /a/2- The assignment of operation 
to input state is random with the constraint that all in- 
put states are used an equal number of times. The ap- 
plication of an operation to its input state produces an 
output density matrix, which we reconstruct using quan- 
tum state tomography [9J [24]. This procedure involves 
nine analysis settings. For each setting, we run the ex- 
perimental sequence 100 times for a total of 900 runs 
per unitary transformation. A single run takes approx- 
imately 37 ms. We compare the measured output state 



to the ideal result using the fidelity [2"5], /(/Oidoai, Pcxp) = 

[Tr ( a/ \/PideaiPexp y/pideai ) ] 2 ■ Figure 2 shows four exam- 
ples of the output states, and Figure 3a gives a histogram 
of the 160 state fidelities. The 160 operations have a 
mean state fidelity of (/) = 79.1(4.5) %, where the error 
bar is the standard deviation of the measurements. Nu- 
merical estimates [TO] indicate that 3.4% of this distribu- 
tion arises from statistical fluctuations in photon counts 
used in state measurement. We attribute the remaining 
distribution to variability in each operation's suscepti- 
bility to experimental noise. We observe no correlation 
between output state fidelity and input state, as demon- 
strated in Figure 3b, c. The mean output state fideli- 
ties from operating on the 16 input states are distributed 
with a standard deviation of 1.5%, as we would expect 
for the means of 10 measurements which themselves have 
a standard deviation of 4.5%. The primary fidelity loss 
mechanisms are percent-level intensity fluctuations in the 
Raman light ficlds[10 and spontaneous emission |26j; the 
fidelities observed here are consistent with those demon- 
strated previously |10j after accounting for the increased 
number of gates. 

As a further check on input-state independence, we 
conduct quantum process tomography [51 [TO] on 11 of the 
operations. Process tomography reconstructs the com- 
pletely positive linear map £ that describes the qubit 
evolution from initial to final density matrices, pa n ai = 
£(Pinitiai)- The map includes the possibility of exper- 
imental imperfections such as coupling to the environ- 
ment, which leads to nonunitary evolution. We represent 
the map by a 16 x 16 matrix[9] E = J2i j K) 01)- 
For each of the 11 operations, we determine an experi- 
mental process matrix i? oxp and compare it with the ideal 
case by calculating both the entanglement fidelity [27 
F = Tr (-Eideai-Eexp) /16 and the mean state fidelity / 
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FIG. 3: Characterization of arbitrary control, a, His- 
togram of the output state fidelities / for 160 arbitrary uni- 
tary operations applied to one of 16 input product states. The 
mean output state fidelities (/} for each input state b, plot- 
ted showing the standard error of the mean and c, binned 
into a histogram, d, The mean state fidelities / for the 11 
operations reconstructed using quantum process tomography. 
The curves are normal distributions with the same respective 
mean and standard deviation as the data. The narrower dis- 
tributions in histograms c,d are expected because they are 
distributions of means. 



obtained by averaging the output state fidelities from nu- 
merically applying E cxp and -Eidcai to an unbiased set of 
36 input states (formed from the eigenstates of tensor 
products of the Pauli matrices) . These fidelities are re- 
lated by / = (Fd+1) I (d+1), where d is the Hilbert-space 
dimension[27 (here d = 4). The 11 processes have mean 
fidelities of (F) = 73(3) % and (/) = 79(2) %. 

In conclusion, we have demonstrated a programmable 
quantum processor capable of implementing all possi- 
ble unitary operations on two qubits. To address large- 
scale problems, many more qubits and gates will be re- 
quired. In anticipation of such applications, this im- 
plementation used only scalable techniques |10j such as 
long-lived qubit storage, quantum information transport, 
and sympathetic cooling. When implementing a larger 
system, the compound errors from successive operations 
will need to be reduced via error correction [3J. This 
will require much higher gate fidelities than shown here, 
both to achieve fault-tolerance and to reduce error cor- 
rection's computational overhead [28 . Nevertheless, the 
type of device described here could form a processing unit 
in a larger system [5] with programmable registers con- 
nected by multidimensional trap arrays [2T)l I29j or pho- 
tonic networks 1 30 . 



Methods 
Algorithmic details. 

For a given two-qubit operation U, calculating the 15 
single-qubit-gate parameters used in the circuit of Fig- 
ure lb is facilitated by working in the so-called "magic" 
basis [11, 31J given in the main text. Transforming to 
the magic basis from the two-qubit computational basis 
1 1 T) , 1 10) , 1 01} , 1 00) is accomplished by use of the unitary 
matrix 
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(6) 



The properties of the magic basis rely on unit matrix de- 
terminants; thus we first strip U of any global phase by 
dividing it by a fourth-root of its determinant, making 
it a member of SU(4). Global phases exactly vanish in 
any observable quantity, allowing this modification. In 
what follows, matrices in the computational basis are de- 
noted with capital letters, and those in the magic basis 
by lower-case letters; e.g., m — A^MA. 

We first find the three degrees of freedom a, f3, 8 that 
produce the correct local equivalence class. We decom- 
pose the circuit in Figure lb as 

U =(C®D)-V -(A(g>B), (7) 

where V determines the local equivalence class and can 
be generated using the gate operations appearing in Fig- 
ure lb within the dashed box. In order to construct V, we 
transform both V and U into the magic basis as v and 
u and choose a, /3, 8 such that the eigenvalues of vv T 
match those of uu T . (We include a global phase e -i7r / 4 
in V to make it an element of SU(4).) This is done by 
comparing the analytical form of the eigenvalues of vv T 
to those of uu T . Since uu T is unitary, it has complex 
eigenvalues of modulus one: Xj = e l ^ j (j £ {1,2,3,4}). 
We find that a, (3, and 8 are given by the means of pairs 
of eigenvalue phases. One possibility is a = {<pi + </>2)/2, 
j3 = (0i + 4>s)/2, and 8 = (02 + </>3)/2. Since no ordering 
of the eigenvalues is required, there are many such com- 
binations that produce members of C/'s local equivalence 
class. The proof of this assignment is by explicit calcu- 
lation of the eigenvalues of vv T and is analogous to that 
given in ref. [H] for the controlled-NOT (CNOT) gate 
rather than our phase gate. 

Second, we find the four single-qubit rotations 
A, B, C, D that comprise the remaining 12 degrees of 
freedom. Note that w T and uu T are unitary sym- 
metric matrices and therefore have real, orthonormal 
eigenvectors [TTJ [TS]. Because they share eigenvalues, it 
is possible to simultaneously diagonalizc them with ma- 
trices k and / such that 



kvv k 



T7„T 
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(8) 
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Here, k and I are eigenvector matrices whose columns 
have been permuted such that equation ([8| is valid. They 
are both members of SO (4) (if necessary, one of the eigen- 
vectors can be negated to change the matrix determinant 
from -1 to 1). By rearranging equation (|8|, we obtain 

/ = v^k T luu T l T kv* = v^k T lu (v^k T lu) T (9) 

(/ is the identity matrix), from which we define m = 
v^k T lu, also in SO(4). We thus have that 



= l T kvm 



(10) 



where l T k and m are both real and in SO(4). Since they 
are real orthogonal matrices in the magic basis, they rep- 



resent single-qubit rotations. We transform equation ( 10 ) 



into the computational basis and compare it with equa- 
tion Q to find 

5 L> = A (Z T /c) A f (11) 
A® B — AmA^ . (12) 



C 



To finish, we split A ® B and C <g> D into A, B,C,D S 
SU(2) and solve for 0, cj>, and (f> z for each. 
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Supplementary Information 



Example decomposition 



As explicit examples of the decomposition described in the main text and Methods, Supplementary Table 1 gives 
two of the many possible decompositions of the unitary matrix 



U = 



/-0.342 + 0.260i -0.042 - 0.53H -0.690 + 0.125* 0.196 - 0.034i \ 

0.195 + 0.743* -0.138 + 0.090i -0.064 - 0.520i -0.304 + 0.127* 

-0.416 + 0.007i -0.244 + 0.786i -0.332 + 0.012* 0.179 + 0.090* 

\ 0.074 + 0.217* -0.084 - 0.076i 0.320 - 0.145* 0.871 + 0.229*/ 



The process matrix from our implementation of the second of these is shown in Supplementary Figure 1. 
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Supplementary Table 1: Two decompositions of the matrix U shown above. The parameters refer to those shown in Figure lb 
of the main text. The last column indicates an overall multiplicative phase factor to make the composite transformation 
mathematically equal to U, but it has no physical meaning. 




Supplementary Figure 1: Comparison of a, experimental and b, ideal process matrices of U; the fidelities are F = 0.77 and 
/ = 0.82. The process matrix E contains the maps £(\i) for all 16 elements \i) {j\. Since each map yields a 4 x 4 matrix 
£k,i, E is a 16 x 16 matrix labeled by m = 4(i — 1) + k, n = 4(j — 1) + /. For example, the 1 00) (11| (* = 4, j = 1) element of an 
input density matrix is mapped to £ (|00) (llj), a 4 x 4 block of E given by m £ [13, 16] and n £ [1,4]. 
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Details for figure 2 



Figure 2 in the main text shows four examples of the output states formed while characterising the quantum 
processor. Supplementary Tables 2-3 give the input product states and 15 control parameters used to generate them. 
The unitary transformations decomposed in Supplementary Table 3 and used to create the states in Figure 2 are as 
follows. 



/-0.362 - 0.428* -0.501 - 0.125* 

0.396 + 0.534* -0.651 - 0.325* 

-0.255 + 0.350* -0.224 + 0.310z 

\-0.002 + 0.238* 0.204 + 0.130z 



0.054 + 0.083* -0.464 - 0.440*\ 

-0.112 - 0.083? -0.088 + 0.051* 

0.320 + 0.328i 0.529 -0.421* 

-0.744 + 0.456* -0.208 - 0.285*/ 



/ 0.138 + 0.564* 
-0.196 + 0.035* 
-0.054- 0.147* 

\-0.152 + 0.758* 



-0.385 - 0.300* 
-0.256 - 0.388i 
-0.146 - 0.622* 
0.324 + 0.180i 



-0.055 + 0.217* 
-0.576 - 0.027* 
0.706 + 0.113z 
0.277 -0.171* 



-0.611 + 0.004A 
0.386 + 0.512* 
0.226 + 0.075* 
0.289 + 0.273* / 



/-0.645 + 0.445* 
-0.528 + 0.217* 
0.012 + 0.136* 
\-0.079 + 0.187* 



-0.154 + 0.006* 
-0.145-0.151* 
0.534 - 0.495* 
0.602 - 0.201* 



-0.174 + 0.561* 
0.404 - 0.676* 

-0.183-0.016* 
0.037 + 0.031* 



0.127 + 0.037* \ 
-0.050 - 0.092* 

0.250 - 0.596* 
-0.545 + 0.507*/ 



/-0.162 + 0.425* 
0.091 - 0.304* 
0.256 - 0.229* 

\-0.607 + 0.454* 



-0.043 + 0.148* 
-0.614-0.158* 
0.105 + 0.743* 
-0.079 + 0.072* 



0.427 + 0.024* 
-0.372 - 0.433* 

0.209 - 0.369* 
-0.116 - 0.546* 



0.031 - 0.765*\ 
0.104- 0.401* 
0.365 + 0.076* 
0.035 + 0.319*/ 



subfigure 




qubit 1 qubit 2 


a 


(|1> 


-i\0))/V2 |0) 


b 


(|1> 


-i|o»/V5 |o) 


c 




11} 11) 


d 


(|1> 


-i\0))/V2 (|l> + |0»A/2 



Supplementary Table 2: Input states used to produce the states in Figure 2 of the main text. 



subfigure a fi 5 6 A <p A <p zA 6 B <p B Z B 6 C <pc <Pz,c S D <p D <p z D 

phase 

I 2.286 0.840 6.094 1.648 1.402 2.443 1.277 2.627 2.930 5.305 2.593 3.873 4.501 2.938 3.643 i 
b 6.196 1.439 1.119 2.542 3.079 4.151 1.982 1.744 6.140 0.833 2.349 1.947 3.742 0.011 1.487 -i 
c 1.589 5.129 1.721 1.721 1.714 3.344 4.911 0.249 1.107 0.470 0.918 0.339 1.632 0.545 1.607 i 
d 0.389 0.750 1.392 5.423 1.300 6.004 4.824 1.961 3.199 4.351 1.100 1.501 2.352 2.022 1.426 -i 



Supplementary Table 3: Single-qubit-gate parameters used for programming the quantum processor to produce the states in 
Figure 2 of the main text. 



